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Abstract 

We consider a set of N linearly coupled harmonic oscillators and show that the diagonalization 
of this problem can be put in geometrical terms. The matrix techniques developed here allowed 
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for solutions in both the classical and quantum regimes. 
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I. INTRODUCTION 



Quantum mechanics is a theory constructed to deal with isolated systems. However, it is 
impossible to completely isolate a system from its surroundings, and usually it relax towards 
thermal equilibrium with the environment, a phenomenon called dissipation. Otherwise, the 
environment acts over the system due to its thermal fluctuations, leading to the destruction 
of pure quantum states in a much shorter time scale; this phenomenon is called decoherence 

m. 

Nowadays, the main theoretical proposal for describing dissipation and decoherence phe- 
nomena is to consider the system of interest plus the environment as a closed system and 
to determine the effective dynamics of the system of interest. This may be done exactly by 
projection technics ^J, and the resulting equation is very useful for an abstract knowledge 
about the structure of the problem. However, since the exact equation that governs the 
dynamics of a quantum system is non Markovian and non linear, to obtain its solution is 
possible only in some very particular cases. 

A model which has been very explored in the literature is the linearly coupled oscillators 
model proposed by A. Caldeira and A. J. Legget 0, Q|. In this model, the system of 
interest is a harmonic oscillator and the environment is modeled as huge set of oscillators in 
thermal equilibrium in a given temperature; the system-environment interaction is given by 
a linear coupling in the position of the oscillator of interest. The Caldeira-Legget model has 
been very successful in various areas of physics, particularly in quantum optics. As is well 
known, despite its quantitative success in the analysis of several experiments [?], 0| , a sound 
derivation and understanding of the hypothesis involved in its derivation remain to this day 
an open question. 

The objective of the present contribution is to understand in detail the classical and 
quantum symplectic structures of a system composed by an oscillator linearly coupled to 
N — 1 others without tracing out the N — 1 states which represent the environment. In par- 
ticular we show that the exact problem possesses a simple physically transparent geometrical 
structure. In Section 2, we treat the classical oscillators. The N = 2 case is analyzed first; 
this problem, although very simple, is structurally rich, and clearly exhibits the geometrical 
nature of the solution: we analyze sections in the phase space which guide us in a sequence 
of variables transformations (a squeezing, a rotation and another squeezing) that give us the 
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system's normal modes. These transformations encountered for N = 2 are then used to find 
analogous transformations for the generic iV case. In Section 3, we work with the quantum 
oscillators. The transformations encountered for the classical case are used there in order to 
find and analyze the normal modes of the system. Throughout this work we use matricial 
technics which have been very useful to express the algebra involved here in a clear and as 
short as possible form. In Section 4, we comment on our results and their relations with 
modern problems. 



II. CLASSICAL SYSTEMS 



A. Two Coupled Oscillators 

Let us consider the two coupled oscillators' Hamiltonian 

where qi and Pi are canonical conjugate variables, g and Ui are real numbers, g < and 
uji > 0. This may be concerned with the system represented in Fig. for example, where 
we have two oscillators with masses rrii and elastic constants ki coupled (position-position 
coupling,) by a spring whose elastic constant is k. If we define as the displacement of mass 
rrii from the stable equilibrium position, and as its conjugate momentum, the Hamiltonian 
for this system will be 

tj Pi P2 , fo. 2 , ^2 2 , & / \2 

n = Qi H Q9 H — Q2 — Qi) , 

2m x 2m 2 2 Ml 2 ^ 2 2 VH H J ' 

which may be put in the form of equation (P) by setting 



/ / 7 7\\— / ^ ^ 

ai = (rrii + fc)) 4 , u)i = \ , g = . 

q { = a^i and Pi = a~ % Pi. 
The Hamiltonian (JJJ may be written as a product of matrices: 

H = x T • H • x, 
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where 



H 



/ ui g \ 

g u 2 

wi 

\ w 2 / 



<?1 92 Pi P2 

Consider a variables transformation determined by the matrix M (whose inverse exist). In 



order to preserve Hamilton's equations, M must obey the symplectic condition 



3: 



M • J • M T = J, 



where 



J = 




is a N x N null matrix and I is a N x N identity matrix (here N = 2, in Section 2.2 N 
will be generic). In this new variables, the Hamiltonian is given by 

H = x T ■ H ■ x, 



where 



M x and He (M t ) 1 H M" 1 



H 



(2) 



It will be very easy to find the equations of motion for the system if M is such that H has 
the form 

/ tt+ \ 


tt + 
\ 

We will search this M by performing sections in the phase space that shall guide us in a 
sequence of canonical transformations. The matrix M is concerned with the whole trans- 
formations carried out in the sequence, and it will be the product of the matrices of each 
individual transformation. Since the following individual transformations obey the symplec- 
tic condition, M obeys too. 



4 



1. The Transformations 



a. SQUEEZING IN p\ x p 2 : Fig. © exhibits the shapes of two sections in the 
phase space. 

The ellipse in the q\ x q 2 plot is inclined due to the term gqiq 2 , which could be eliminated 
by a variables transformation equivalent to a rotation in the plane q\ x q 2 . However, the 
symplectic condition obliges us to make the rotation in the plane pi x p 2 too, and it would 
produce a term proportional to pip 2 . The solution is to perform a squeezing in the plane 
Pi x p 2 first, turning the ellipse into a circle. Thus the rotation we shall make later will not 
cause a new coupling term. 

The matrix that carry out this squeezing is 

fa 

a' 1 



Mi 



\ 




a" 1 
\ a ) 



where 



Defining 



we may write 



a = 



u 2 




y / LUiLU 2 , &i = \ — and Q 2 = \ — 




H^Mf^-H-Mr 1 





( fii 
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1 
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n 2 








2 
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The new variables are given by 



Pi" 



.(1) 



Mi.x, 



5 



and the Hamiltonian may be represented in the form 

H= (x«) T -H«-x« 



The sections will be like the ones in Fig. (J3J). 

b. ROTATION IN q± x qlp : Now the rotation will not create new nondiagonal 
terms. The transformation matrix is 

' cos^ simp ^ 

- sin ip cos if 

cos <p sin ip 

\ — sin <p cos if J 

and we define 



M 2 = 



where 



H( 2 ) = (M^-H^-M- 1 



1 f Hi? H§ 

2 I Hg> Hg 



H 



(2) 
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Qi cos 2 + Q 2 sin 2 ip + 2g sin ip cos 93 sin <p cos (/? — &i) + 9 (cos 2 — sin 2 ip) 
sin </? cos <p (O2 ~~ ^1) + <7 (cos 2 </? — sin 2 (pj fli sin 2 93 + Q 2 cos 2 tp — 2g sin y9 cos </> 



We now choose 



H (2) _ jj(2) 
X1 12 ~~ ""-21 



001 ™ / oj 

tt(2) 
) X1 22 







uj 



cos (p 



sin </? 



\ 



2 |1 + 



(Q 1 -Q 2 ) 2 + 4g^ 



V 

in order to make diagonal. Thus 



H< 2 > = - 
2 



/ cj+ \ 

w_ 

w 

\ u ) 



where 



10. 



2 ( («i + ^ 



Defining the new variables as 



(fli-n2) 2 + V , 



(O x -0 2 ) 2 +4^ 



ri 2) 

U 2) / 



X ( 2 ) = M 2 -x«, 



we may write the Hamiltonian in the form 

For these rotated variables, the sections will be like the ones in Fig 

c. SQUEEZING IN qf ] x p[ 2) AiVX) /AT g^ 2) x pf : The matrix H^ 2 ) is not in 

the final form (equation (J2J)), because the coefficient of each coordinate is not equal to the 

coefficient of its conjugate momentum, as may be seen in Fig 
The following transformation does the adjustment: 

/ a + \ 



M, = 



a- 
aT 1 



\ aZ 1 J 



where 



-(^) 

We find the matrix H (equation (J2J) by defining 



1/4 



H = (Mj ) 1 • H^-Mg 1 

/ n+ \ 

1 

2 Q + 
\ / 
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where 



The matrix which corresponds to the set of transformations carried out is 



M = M 3 ■ M 2 ■ M 



i- 



Thus we may write 



H= (M T ) 1 h m 1 , 



EE X EE M X. 



and the new variables are 

' Qi 
Pi 

W 

The expression for the Hamiltonian using x is 

H = x T • H • x. 

Fig. (P) shows the equality of the coefficient of each coordinate and its conjugate mo- 
mentum. 

The parameters Q± are given in terms of the original parameters by 



The constant Q + is always real, since ia% and u>2 are positive real numbers. The constant Jl_ 
will be real when 

wiw 2 > g 2 - 

If we are dealing with a spring- mass system like the one in Fig. (0), this condition is always 
satisfied, since it corresponds to 



(fci + k) (k 2 + k) > k 2 
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2. The Equations of Motion: 



Since M refers to a canonical transformation, the Hamilton's equations are preserved. 
Thus we may write them, in the matricial form M, for the variables in x: 



^ j -J^( j '!)'(* r ' fl '*) 



where 



It yields 



d_ 
<9x 



JL Odd 
dq\ dq2 dpi dp2 



( fax \ ( \ 



I?2 

V Ip2 / 



whose solution is 
/&(*)\ 



\P2(t)J 



V 



Pi(0) 

Wo)/ 



/ cos(fi+t) sin(fi + t) \ 

cos(fi_t) sin(fi_t) 

-sin(fi+t) cos(fi+t) 

-sin(Q_t) cos(ft_t) / 
It is clear in the equation above that the frequencies of the normal modes are fl + and fi_. 
Defining 

/ cos(fi + t) sin(fi + t) \ 

cos(fi_t) sin(fi_t) 

sin(fi + t) cos(fi+t) 
-sin(n_i) cos(fi_t) / 
we may write the equations of motion for the original variables: 



N 



V 



x(t) = M _1 -x(t) = M _1 -N ■ x(0) = M" 1 ■ N ■ M ■ x(0). 



B. N Coupled Oscillators: 

Consider now one oscillator linearly coupled to iV — 1 other oscillators. Of course the 
system we just solved is a particular case of this system. Thus the transformations we found 
there are useful to suggest the form of the transformations to be effected here. 
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The N coupled oscillators' Hamiltonian is 

j N N 

i=l i=2 

where gj and pj are canonical conjugate variables, all the Uij are real numbers, u;^ > and 
Ui j ^ for % ^ j. This Hamiltonian may be written in a matricial form: 



H 



1 




where 



Hq = 



^1,1 ^1,2 " ' " ^1,AT- 
^1,2 ^2,2 

; o ••. 







^l,N-i '■ 

u hN 



wn-i,n-i 







^N,N j 

\ 








uj lt i • • 

U 2y 2 

; o '•. '•. ; 

I •• LUn-1,N-1 

• • • U>N,N 

q 1 q 2 ■■■ q N -i Qn 

Pi P2 ••• Pn 



i Pn ) • 



The aim of the following transformations is to diagonalize the matrix 

Hq 

Hp 

All of them are canonical, as may be verified using the symplectic condition, and reduce to 
the transformations found in Section 2.1 when N = 2. 



1. The Transformations 



a. SQUEEZINGS INP: Regard 

M s = Si • S 2 • ... • Sjv„i • Sjv, 
10 



where 







Si = (w M ) 4 





; o 



\ 



; ••. o : 
o ••• o '•. o 
oo ••• 1 



= (^,0 4 (^hi) 2 > [Si\j,k = (^,i) 4 $j,k for 3 7^ Tne variables transformation to be 
effected is 

' Q s \ = / M s W Q 
Ps )~ { M^ 1 J \P 

and the Hamiltonian may be written as 



ft i)( 



Qs 

Ps 



where 



H Qs 
H Ps 



M s 
M" 1 

Mg 1 • H Q M 1 




H Q 

Hp 

s 

M s Hp M s 



The matrices Hq s and H Pg are given by 

/ 



M s 
Ms 1 



H Qs 



II 



Ps 



01.1 01,2 ' ' ' 9l,N-l 9l,N 

01.2 02,2 ••• 

: ■•• ••• : 

' • 9n-i,n-i 



\ 



9l,N-l '■ 

\ 9i,n ••• 

/ G ••• 

o g '•■ ; o 

: ■•■ ■•■ : 

••• G 

y ... G 







9n,n J 



(3) 
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with 



N 



G = n, 



1,1 1 



1=1 



9ij 



Uj,j^LOi,iUJj,j 

G 



This transformation corresponds to N squeezings, each one involving and Pi. Notice 
that G and g^j are real numbers (G > 0, > and g^j ^ for i ^ j). 

b. ROTATIONS IN Qs- - Since Hq s is a real symmetric matrix, it has real eigen- 
values and may be diagonalized by a real orthogonal matrix M R . Thus we define 



Hq r = M 



R ' H Qs ' M R 



(4) 



Ai,i ••• 
\ 2 , 2 



V 

where the A,- , are real numbers and 








AjV-l,iV-l 

A 








N,N 



Mr • M R = I, 



(5) 



(I is the N x N identity matrix). To find the \ iti may be non trivial, since, to this end, it 
would be necessary to solve a N degree polynomial. We do not have a general expression 
for M R , but we show in the Appendix that it is possible to write M R as the product 



M R — Rn 1 ■ R-N-2 • ••• • R2 • R-l, 



(6) 



where 



R; = R^jv ' Rj,W-l ' ••• ' Ri,i+2 ' Rj,i+1, 

and the Rjj are matrices whose elements [Ri,j] feJ are given by 











i,j 


= sinctij, 


[Ri,j 


'3,i 


= — sin aijj , 






= 1 for k 7^ % and k 7^ j, 




k,l 


= for the other cases. 



12 



We may interpret TU,j as rotation in a N dimensional space, and Mr as a sequence of such 
rotations. 

An intuitive view of the Mr decomposition may be achieved. Equation (JHJ) permit us to 
interpret the Mr lines as the coordinates of vectors which form an orthonormal basis in a 
TV dimensional space. Of course this coordinates are written in another orthonormal basis, 
which we call the canonical basis. If we apply a suitable sequence of rotations (represented 
by M R ) in the basis represented by Mr, we can make it coincident to the canonical basis, 
and it will be represented by an identity matrix (M R • Mr = I). In the N = 2 and N = 3 
cases it is not difficult to imagine the rotations being performed. 

The transformation to be effected here is 



Ms 1 

We may write the Hamiltonian as 






2 lQ*p R j^ Hpr 



> 


II 


' Qr \ 


H PR 1 




V PrJ 



where 



H 



Qr 







H Pr 




c. SQUEEZINGS IN Q R x P R . Consider now 

M T = Ti ■ T 2 • ... ■ Tn-i ■ Tn, 



where 



/ 



Ti = (a> M )' 



10 ••• 

o '•. '•. ; o 

■■• '•• 

••• 1 



\ 
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Defining 




M T 
M T X 



Mrp 1 • H Qr • M T ! 









M T ! 



we have 



/ 



H 



Qt 



H 



p T 



o n 



2,2 




















ooo n N - 1)N -! 

\ 



\ 





Qn,n J 



where 



Of course M T corresponds to N squeezings like the ones carried out in Section 



2.1. 



The final variables will be given by 



M T 
M T X 



and we may write the Hamiltonian as 









Mg 1 




^(^)( H QT H° P 




2. The Equations of Motion: 

Since the system is decoupled for Q T and P T , the equations of motion are easily calcu- 
lated: 

' Qt(*) \( Ac A s \ f Q T (0) 
P T (t) ) \ -A s A c ) \ P T (0) 

where A c and A s are N x N diagonal matrices whose elements are 

[ A c] M = cos(a,^), 
[A s ] = sin{Q hl t). 
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The frequencies of the normal modes are the Vt^. 
Notice that 

I Ac A s \ 
A = \ = Ai • A 2 • . . . • A N _x • A N , 

\ -A s A c J 

where Ai are 2N x 2N matrices whose elements are given by 

[Ai] M = [Ai\ N+i:N+i = cos(ft M t), 

l A i\i,N+i = sin(fi M *), 

[M N +i,i = -ffln(fi M t), 

[Ai] fe k = 1 for A; 7^ i and k ^ N + i, 

[Ai] feJ = for the other cases. 

In this manner we may see A as a sequence of iV rotations in the planes qi T x p iT , where 
q iT and p iT are the elements of Q T and.P T : 



Qt — ( Qi T Q2 T ■ ■ ■ Qn-i t Qn t ) ) 
p t = (pi T P2 T ■■■ Pn-i t Pn t ) • 



Defining 





M T \ / M R \ / M s 
M T x y ^ Mr y y M s ! 

the equations of motion for the original variables may be written as 



Z 1 A Z 



III. QUANTUM SYSTEMS 

As is well known quantum problems involving harmonic oscillators and linear couplings 
are in general equivalent to the classical ones, except for essentially quantum initial condi- 
tions. As we show here the transformations are the same as the ones obtained in the classical 
case. 
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Consider a quantum oscillator linearly coupled to N—l other oscillators. The Hamiltonian 
operator in the so called rotating wave approximation (RWA) shall be written as 

N N 

H = ^g iti a\ai + ^T,9i,i (a\a,i + 0J01J , (7) 

i=l i=2 

where > 0, g it j ^ for i ^ j, and we set H = 1. The operators and a\ are bosonic 
creation and annihilation operators: 



Si,j. 



(8) 



This Hamiltonian may concern a microwave mode constructed in a superconducting cavity 
which is connected to iV — 1 other superconducting cavities by waveguides, as it was assumed 



in Ref. [lCj for N = 2, where no environment was considered. Otherwise, we may regard 
the oscillator related to a\ as the system of interest and may use the other oscillators for 
simulating the environment. 

In the following, the matrices of Section 2 will help us to get information about this 
system. The orthonormality of M R , Eq. (JSJ), will be used in several steps below. 

Defining 



M„ = 



a 2 



/„t\ 



and M a t = 



\ a N J 



the Hamiltonian (|ZJ) may be written as 



H = (M a t) T ■ HQ S -M a , 
where Hq s is given in Eq. (J3j). Using Eq. we write the Hamiltonian in a diagonal form: 



N 



H = (M st ) T ■ H Qr -M, = ^Ai.iotai 



where 



M fl = 



a 2 



y a N J 



( -J \ 



= M R -M a and M at = 



\ a N J 



M R -M 



(9) 
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The operators di and a\ above are the bosonic operators related to the normal modes 
of this system, whose frequencies are A^. They obey the usual commutation relations for 
bosons, as may be seen by first writing the commutation relations for the a, and a] in a 
matricial form, 

M a .(M at ) T - (M at -(M a ) T ) T = I, 

and then observing that 

M a - (M at f - (M at - (M,) T ) T = M R - (m u - (M at ) T - (m oV (M a ) T )^) M R = I. 
If |0) is the vacuum state related to the original operators a iy we must have 

a, |0> = |0> 

for any aj. It is easy to see that |0) is the vacuum state related to the operators aj, since 
the (ii are linear combinations of the a^. 

The number of excitations of the normal modes is a conserved, since 

" N 



i=i 



0. 



Observing that 



N 



N 



= (M a t) T -M a = (M at f -M a = J2 C 



we see that the number of excitations is conserved for the original modes too. 

Let us define 1 1 j) as the state where there is one excitation in the i th oscillator and vacuum 
in the other oscillators. Thus, 

\h) 



\\1n)J 



= M at |0) 



Analogously, we may define one excitation states for the normal modes: 



Mil) 



|i 2 > 

{\In)J 
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= M Bt |0) 



Observing that 



it is easy to see that 



where 



e iHt \li) = e iA ^|li), 
M,!) = M R -M, 1): 



e iHt M|i) = M^-e mQ R t -M R -M|i) 



e iH QR* 



••• 



e iA2 - 2 * 



V 



o 



o ••• e iXN ' Nt 
what gives the temporal evolutions of the states 



IV. CONCLUSION 



We performed a detailed analysis of a system composed by one oscillator linearly coupled 
to N — 1 other oscillators. An intuitive view about this problem may be achieved: for N = 2 
and N = 3, the variable transformations carried out in order to find the system's normal 
modes may be geometrically visualized, and for other values of N analogous transformations 
occur. If we consider that the central oscillator is the system of interest and the other 
oscillators are the environment, we will be able to study decoherence, e. g., comparing 
results obtained for this system to results obtained using master equations |7], or even to 
results coming from the laboratory 

The matricial technics developed here are easier to deal with in the classical systems. 
The classical results can be almost directly used for the quantum systems. Of course, other 
topologies for the oscillators may be investigated using similar matricial technics, e. g., 
two oscillators coupled to iV other oscillators. If we consider the two central oscillators as 
the system of interest and the N other as the environment, conditions for the appearing 
of decoherence-free subspaces |3, ll jl ll^ . jl^ may be searched (two oscillators subjected 
to the same bath were studied in |8j, [lfj [l7( using master equations). Knowledge about 
this topic is important in the quantum information context, since the main problem in 
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this area nowadays is the deleterious action of the environment over quantum coherences 
(decoherence). Detailed comparison from our results with master equation results will be 
the subject of a forthcoming publication. 
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APPENDIX 

In order to see why the decomposition of Section 2.2.1 is possible, we first notice that it 
is possible for N = 2: 



\ — sin 0^2 cos 0:1,2 / 

Then we prove the general case showing that if the decomposition is possible for a (N — 1) 
x (N — 1) orthogonal matrix it will be possible for a N x N orthogonal matrix. 
Consider a N x N real matrix A. Observing that for k 7^ 1 and k 7^ j 




for N = 2. 



[ A - R l-]i, fe = [ A l 



l,k ' 



and that it is always possible to find a value for a±j which yields 




we see that it is possible to chose all the aij in such a way that 



B EE Mr • RJ 2 • ... 



I 1 "O 



1,JV-1 




1,JV 







\ 



V 



Bn,i B Nj2 




We may write 



M R = B • Ri. 



(10) 
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Observing that if A is orthogonal then A • Rf • is orthogonal too, we see that B has the 



B = 



form 

f 1 0---0 


; c 
o 

where C is a (N — 1) x (N — 1) orthogonal matrix. Since we are assuming that the decom- 
position proposed is valid for a (N — 1) x (N — 1) orthogonal matrix, we may write 



V 



C — Dn-2 • Dn 3 



D 2 D l5 



where the Di are the (N — 1) x (N — 1) matrices analog to the Ri. Thus 

d o...o\ I 



: Dn-2 



B 



V 



o 



/ 



1 0---0 


D N _ 3 



V 



o 



Notice that we may regard 



/ 



/ 



Ri+l 



V 



V 

1 0---0 

o 

; Di 

o 





; d 2 





/ 



and then 



B = Rn i • Rn-2 ■ ••• • R3 • R2- 

Using equation (|10jl we have Mr in the form given by (JHJ). 



' 1 0---0^ 
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FIG. 2: Sections in the phase space: q\ x q2 and p\ x p2- 
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FIG. 6: Sections in the phase space: gi x pi and (72 x j?2- 
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